C 13, </</ /<£T""~ „ aja i 
NBS MONOGRAPH 85 



A Study of Lunar Surface 
Radio Communication 




Documents Collection 
OCT 19 1964 



U.S. DEPARTMENT OF COMMERCE 

NATIONAL BUREAU OF STANDAROS 



THE NATIONAL BUREAU OF STANDARDS 



The National Bureau of Standards is a, principal focal point in the Federal Government for assuring 
maximum application of the physical and engineering sciences to the advancement of technology in 
industry and commerce. Its responsibilities include development and maintenance of the national stand- 
ards of measurement, and the provisions of means for making measurements consistent with those 
standards; determination of physical constants and properties of materials; development of methods 
for testing materials, mechanisms, and structures, and making such tests as may be necessary, particu- 
larly for government agencies; cooperation in the establishment of standard practices for incorpora- 
tion in codes and specifications; advisory service to government agencies on scientific and technical 
problems; invention and development of devices to serve special needs of the Government; assistance 
to industry, business, and consumers in the development and acceptance of commercial standards and 
simplified trade practice recommendations; administration of programs in cooperation with United 
States business groups and standards organizations for the development of international standards of 
practice; and maintenance of a clearinghouse for the collection and dissemination of scientific, tech- 
nical, and engineering information. The scope of the Bureau's activities is suggested in the following 
listing of its four Institutes and their organizational units. 

Institute for Basic Standards. Electricity. Metrology. Heat. Radiation Physics. Mechanics. Ap- 
plied Mathematics. Atomic Physics. Physical Chemistry. Laboratory Astrophysics.* Radio Stand- 
ards Laboratory: Radio Standards Physics; Radio Standards Engineering.** Office of Standard Ref- 
erence Data. 

Institute for Materials Research. Analytical Chemistry. Polymers. Metallurgy. Inorganic Mate- 
rials. Reactor Radiations. Cryogenics.** Office of Standard Reference Materials. 

Central Radio Propagation Laboratory.** Ionosphere Research and Propagation. Troposphere 
and Space Telecommunications. Radio Systems. Upper Atmosphere and Space Physics. 

Institute for Applied Technology. Textiles and Apparel Technology Center. Building Research. 
Industrial Equipment. Information Technology. Performance Test Development. Instrumentation. 
Transport Systems. Office of Technical Services. Office of Weights and Measures. Office of Engineer- 
ing Standards. Office of Industrial Services. 



* NBS Group, Joint Institute for Laboratory Astrophysics at the University of Colorado. 
** Located at Boulder, Colorado. 



UNITED STATES DEPARTMENT OF COMMERCE «> Luther H. Hodges, Secretary 
NATIONAL BUREAU OF STANDARDS • A. V. Astin, Director 



A Study of Lunar Surface 
Radio Communication 

L. E. Vogler 




National Bureau of Standards Monograph 85 
Issued September 14, 1964 

For sale by the Superintendent of Documents, U.S. Government Printing Office 
Washington, D.C., 20402 - Price 70 cents 



Library of Congress Catalog Card Number: 64-60070 



CONTENTS 

PAGE 



1. 


INTRODUCTION. ....................... 


o 1 


2. 






3. 


GROUND WAVE ATTENUATION A t ......... i ........ . , 


. . 3 


4. 


ANTENNA EFFECTS. ........................... 


. 19 


5. 


NOISE EFFECTS .............................. 


25 


6. 


POSSIBLE IONOSPHERIC EFFECTS „ ................... 


29 


7. 


CALCULATION OF REQUIRED POWER .................. 


36 


8. 




, 48 


9 e 


ACKNOWLEDGMENTS . . . 


48 


10. 


APPENDIX: INPUT IMPEDANCE 


. 49 


He 


REFERENCES 


, 55 


12. 




. 57 


13. 


GRAPHS OF INPUT IMPEDANCE 


. 87 



iii 



A STUDY OF LUNAR SURFACE RADIO COMMUNICATION 



L, E. Vogler 

The problem of point -to -point radio communication on the moon is discussed, 
and equations and curves are presented to estimate power requirements in lunar com- 
munication systems „ Assuming a smooth surface, consideration is given to ground 
wave attenuation over both layered and non-layered grounds, antenna ground losses in 
situations where ground screens are impractical, noise level estimates in the receiv- 
ing system, and the effects on propagation of possible lunar ionospheres. An example 
of the calculation of required power for a particular communication system is given, 
and further studies are suggested. 

1. INTRODUCTION 

This monograph presents the results of a study investigating the problem of point-to-point 
radio communication on the surface of the moon. It is intended as a survey of all the work accom- 
plished in a project sponsored by Jet Propulsion Laboratory and, consequently, much of the material 
has appeared in previous papers concerned with various phases of the study [ Vogler, 1963a, 1963b; 
Vogler and Noble, 1963, 1964] ; the rest of the material, notably the discussions on the effects of 
layered grounds and possible lunar ionospheres , has not been published before. 

No attempt is made here to analyze in detail one particular type of communication system; 
rather, it is hoped that this study may provide part of the general information necessary to the design 
of any future lunar surface communication link. For example, the antenna used in a permanently 
based installation would most likely be much more elaborate than the one needed for communication 
by a mobile exploring party. Once the purpose of the system is determined, the curves herein may 
be used to help estimate the frequency, range, and power requirements. 

Most of the curves, aside from those figures describing illustrative examples, are plotted in 
parametric terms for the s imple reason that very little information exists at the present time con- 
cerning the exact nature of the lunar envir onment; however, as more and more data become available , 
the curves, using the better estimates of parametric values, will still be useful in the analysis of 
communication requirements. 

The study is primarily concerned with the propagation aspect of lunar surface c o mmuni c a ti on , 
although considerable emphasis has been placed on the evaluation of antenna power losses due to the 
proximity of the ground. For low antenna heights and poorly conducting soils, these losses may be 
far from negligible, especially where ground screens are impractical. The effects of noise and its 
relationship to the desired signal is discussed briefly, consideration being given to galactic and 
receiving system component noise. 

The lunar model assumed is that of a smooth sphere of radius r = 1738 km. The effects of 
rough terrain on s ignal propagation are not included in this analysis, although it is recognized that 
such things as knife-edge diffraction, "obstacle gain", etc. , might be important at some locations 
on the moon's surface; a study including these mechanisms could be undertaken for particular surface 
areas of interest. The present curves are expected to be most useful for applications in the frequency 
range from about 1 kc/s to 10 Mc/s. Throughout the present analysis, antenna heights are assumed 
to be low enough such that height gain effects are negligible. 
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2. SYSTEM LOSS 

The concept of system loss [Norton, 1959; Wait, 1959; CCIR, 1963 ] will be used to 
describe the effects of the various system parameters. System loss L g is defined as the ratio, 
expressed in decibels, of the radio frequency power input to the terminals of the transmitting antenna 
p^_ and the resultant radio frequency signal power available at the terminals of the receiving 

antenna p : 

r 

L s = 10 log(p t /p r ) = P t - P r . (2.1) 

Unless otherwise specified, capital letters will be used to designate the decibel equivalent of the 
corresponding lower case letters; e. g. , = 10 log ^ p^ . 

The system loss may be divided further into a free space component that would be the loss 
expected between isotropic antennas situated in free space, an attenuation relative to a free 

space inverse distance field which accounts for the effects on propagation of the intervening terrain, 
and the transmitting and receiving antenna components L^, and describing the effects of 

the particular antennas used. Thus , (2. 1) may be written as 

P = L + P =10 log(4iTd /\) 2 + A- G+L+L+P , (2. 2) 

tsr to o tptrr 

where d is the arc distance separating the two antennas and X. is the free space wavelength 
o 

expressed in the same units as d^. The term denotes the path antenna gain which is a function 

of the transmitting and receiving antenna maximum directivity gains relative to isotropic antennas , 

G and G , and any aperture -to -medium coupling loss or polarization coupling loss . For the 

t r 

purposes of this analysis , the latter two losses will be assumed negligible such that G^ = G^ + G^ ; 
however, it should be noticed that, in general, the path antenna gain will be less than or equal to the 
sum of the maximum directivity gains. The symbols and L,^ denote antenna losses arising 

from the antenna's proximity to the ground: ^ = 10 log (r /r f ) , where r is the actual input 

resistance and the - radiation r esistance in free space. Notice that, in order to calculate trans - 

mitter power, any transmission line loss and appropriate circuit losses should be added to the right 
side of (2.2). 

The power available at the receiver consists not only of the desired signal but also of radio 
noise arising both from within and from without the receiving system. If r (s /n) denotes the desired 
minimum signal-to-noise power ratio that will provide a given grade of reception as measured at the 
receiver predetection output, and an operating noise factor f^ is defined as the ratio of ( 1) the 
signal to a reference Johns on -noise power ratio that is available from a loss -free receiving antenna 
to (2) the receiver predetection output signal-to-noise ratio r (s /n) , then p_^ may be expressed as 

P r = £h B t o h r ( s/n ) /je c > ( 2 - 3 ) 



^The relationship between f and the effective "noise temperature"' T gff used in some discussions 
is given by T q££ = t Q f. 
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where k t b is the reference Johns on -noise power. The symbol k denotes Boltzmann's 
constant ( = l o 38054 X 10 joules / degree) , t is a reference temperature in degrees Kelvin, 
b is the effective noise bandwidth in cycles per second, and I is the ratio of the antenna input 
resistance to its radiation resistance 

Using the above definitions and assuming I cr» . 



(2 C 2) may now be expressed as 



= 32 c 45 + 20 log d (km) + 20 log f + A - ( G, + G ) + L +R + F + B + 10 log (k t ) , (2 e 4) 
t ° o ° mc t t r t & B o 

where d^ is measured in kilometers, f^^ denotes the radio frequency in megacycles per second, 
and R p 10 log r (s /n) . The following sections will discuss the various components of (2. 4) . 

3o GROUND WAVE ATTENUATION A 

Numerical procedures for the calculation of electromagnetic fields diffracted around a smooth 

homogeneous sphere have been developed by various authors [ Norton, 1941; Bremmer, 1949 ] . 

The ground wave may be expressed as a series of residues that depend on the radius of the sphere 

r , the arc distance d , the antenna heights h, and h , the free space wavelength X, the 
o o to 1 2 r 6 

relative dielectric constant and conductivity cr of the ground, and the polarization of the wave. 

As stated in section 1, the present analysis assumes negligible height gain effects corresponding to 

antenna heights near zero; for elevated antennas, estimates of height gain may be obtained for beyond 

the horizon paths from another paper [ Vogler , 1964 ] „ 

In parametric form the ground wave attenuation relative to a free space inverse distance field 

o , 

A^_ is conveniently plotted as a function of three parameters : K, b , and . Through the recent 
work of Wait [ 1962 a ] both non-layered and layered grounds may be accounted for by generalizations 

o 

of the definitions of K and b . Thus, with the subscripts h and v referring to horizontal and 
vertical polarization respectively, 



K h = 



(2-irr /X) 3 I T, 



K 



^(2irr o /\) 1 



Q 



(3. 1) 



where 



Tj =|"<« r - 1 > 2+l 



% 



(3.2) 



and 



o o -l/ /6 r^ S \ ° o - \ f 

b h = 180 - tan [ ——J - , b^ = 2 tan ( — } - tan 



-If r 



(3. 3) 



with 



-7 4 
2X10 c X cr ^ 1.8 X 10 (r(mhos/m)/f 



3 



Both and s are dimensionless and refer to the top ground layer unless otherwise specified. 

The parameters X and cr are in meters and mhos per meter respectively, and c is the 
velocity of light. 

. . o 

The symbols |Q| and q appearing in (3. 1) and (3. 3) denote the magnitude and phase 
of a "stratification factor", Q == |q| exp (i q /2) , that accounts for the effect of layering in the 
ground. In general Q is a complex quantity dependent on the heights and electromagnetic ground 
constants of the various layers, whereas for a homogeneous (non-layered) ground, Q is equal to 
unity. Wait [ 1962b] has derived formulas for Q for the general case of n layers and has pre- 
sented curves [ Wait, 1962 a] for two and three layer s ituations ; however, unless specific values 
for the heights and ground constants of the various layers are given, a graphical representation 
becomes quite complicated for a multi - lay e red ground. For the present discussion, restriction to 
vertical polarization over a two-layered model will be assumed. 

Letting the subscripts 1 and 2 refer to the upper and lower layers respectively, may 
be expressed as 



/ \ i(b° -b° )/2 ( i(270°- b° )/2. 

(l T vJ/lTjV V1 Vj +tanh{(Z, V X)lT hi |e^ ^ } 

/ x i(b° -b° ]/Z , i( 270°- b° /2. 



where is the height of the upper layer measured in the same units as \ a | T^ | and | T^ | are 

given by (3. 2) us ing the appropriate values of € and cr characterizing layer 1 or layer 2; 

o o r o 

b^ and b^ are obtained from ( 3. 3) after setting q = 0. As 1^ becomes very large (corre- 
sponding to homogene ous ground), approaches unity and | | =1, q^ = . 

Graphs of I T, I and I T I as functions of e and s are shown in figures 1 and 2 
1 h 1 1 v 1 r o 

where the expressions for ^ refer to a homogeneous ground ( |q| = 1) ; similar plots of b^ 

o ' o 

and b^ are given in figures 3 and 4, again for homogeneous ground (q = 0). It can be seen that 
b° is always positive for either polarization as long as no layering occurs. 

o 

The variation of and b^ for a two-layered ground is shown in figures 5 through 8 . 

The curves were calculated from ( 3. 1), ( 3. 3), and ( 3. 4) using the values of the ground constants 

and upper layer heights indicated. Only value s of K for 0.01 < K < 10 are plotted since the 

o V 

variation of as a function of b is negligible outside this range (see figures 9 - 14) The curves 

o o 

of b^ indicate that, as the conductivity of the lower layer cr-, increases, b^ is more likely to 
become negative. Thus, if it is found that the moon consists of a highly conducting core underlying a 
thin poorly conducting crust, ground wave communication could encounter "fade-outs" for certain 
frequencies and ranges. This is indicated in figures 9 through 14 which show the ground wave atten- 
uation for zero antenna heights versus the distance parameter x^ , where 

JL 

x' = f 3 d (km) . (3. 5) 

o mc o 
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€ r 2 +S 2 



1/2 



K v s [(2 7T r /X ) 1/3 1 T v |] " = 3.02 x IO- 2 /|T v | - f mc 
r : MOON'S RADIUS = 1738 km. 
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FOR LARGE S= |T V |~IA/S 



U 0.0I 
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Figure 2. The parameter } T ) for vertical polarization. The 
formula for K assumes a homogeneous ground. 
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Figure 6 e and for two layers Q upper layer height: 10 m 8 
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Figure 7 @ and for two layers 9 upper layer height: 100 m. 
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Figure 10 o Ground wave attenuation for zero antenna heights, 
b = -45° . 
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Figure 13. Ground wave attenuation for zero antenna heights, 
b = 90° . 
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Figure 14® Ground wave attenuation for zero antenna heights, 
b = 180 @ . 
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In figures 9 (b = -70 ) and 10 (b = -45 ) the curves show fading effects over short intervals of 

! 

for certain value s of K D As b becomes positive, these fade-outs disappear and is seen 

to increase smoothly with distance. 

The calculations of the attenuation curves were not done by numerically evaluating the complete 

residue series, since a prohibitive number of terms are necessary for small values of x' , In this 

o 

region (the left hand portion of figures 9 - 14) use was made of the modified flat earth formula 
developed by Wait [ 1962b, see equation ( 108) ] in which the leading term is the plane earth attenua- 
tion function and successive terms account for the curvature of the sphere. The formula with two 
curvature correction terms was employed wherever valid until x^ became too large, at which point 
the first term of the residue series could be used from there on out. 

Equation (2. 4) is strictly applicable only for the case of antennas separated a sufficient 
distance apart (d >> \) such that the induction and static fields may be neglected. Because of 

this restriction, A is not shown in the figures for values of x' < 1 . As d (or x' ) goes to 
t o o o 

zero, A^_ approaches 20 log( 1 /2) = -6.02, which would correspond to the surface wave field intensity 
expected between short vertical electric dipoles situated near each other on a perfectly conducting 
plane. 

It is apparent from figures 9-14 that the attenuation is less for the larger values of K, 

these values corresponding to vertically polarized waves at the lower frequencies. For example, 

for K > 0. 1 and f < 1.0, the attenuation is less than 50 db out to distances of 100 km. On 
mc 

the other hand, with horizontal polarization, is usually much less than unity for all frequencies 

and the attenuation is greater. 

4. ANTENNA EFFECTS 
The terms ^ and = + (the subscripts t, r referring to e ithe r transmitting 

or receiving antenna) in (2. 2) describe the effects of the particular transmitting and receiving 

antennas used in the communication system under consideration. L , called the ground proximity 

t, r 

loss, is defined as the ratio expressed in decibels of the input resistance of the antenna r to its 

free space radiation resistance r^ . As the height of the antenna above the surface is increased, the 

ratio of the resistances approaches unity and I^_ or effectively become zero. For heights 

near the surface, the input resistance is a function of the electromagnetic ground constants e and 

cr i the radio frequency f^^ , and the antenna height h above the surface. 

G^ ^ denotes the free space maximum directivity gain of an optimally oriented loss-free 

antenna above an isotropic antenna, e, g. , G^ = G^ = 10 log (3/2) = 1.76 for elementary dipoles and 

G = G = 2.15 for half -wave dipoles . 
t r 

Because present estimates indicate the conductivity of lunar surface material to be relatively 
poor, a numerical evaluation of input impedance was undertaken [ Vogler and Noble, 1963, 1964 ] 
considering four types of antennas: vertical and horizontal electric and magnetic dipoles. Short 

dipoles of length di are assumed, the current in which varies as I e 1Wt with I a constant 

o o 

along the antenna; the magnetic dipoles are then equivalent to small circular loop antennas of area 
d A, the axis being in the direction of the dipole. The ground beneath the antenna is assumed to be 
homogeneous and characterized by its relative dielectric constant e and conductivity o" <> 
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The input impedance Z, obtained in terms of the field at the antenna, may be divided into a 
free space component and a component AZ describing the effect of the presence of the ground: 



Z = Z f + AZ 



(4. 1) 



Explicit expressions for the impedance change due to ground normalized by the free space radiation 
resistance, A Z/r , where VED, HED, VMD, HMD refer to the initial letters of the four types 
of dipoles considered, are then given by 



VED: AZ/r f = i(3/2 a ) 



HED: AZ/r f = i(3/4 a ) 



VMD: AZ/r f = i(3/2 a ) 



HMD: AZ/r f = i(3/4 a ) 



^(N 2 ) + I Z (N 2 ) 1 , r f = 20 (3 2 (di) 2 , 



1^1) + I 2 (N ) 



+ I 2 (l) 



I^N ) +I 2 (1) 



r f = 20 p 2 (di) 2 , 



, r = 20 (3 4 (dA) 2 , 



, r f = 20 (3 4 (dA) 2 , 



(4. 2a) 
(4. 2b) 

(4,2c) 
(4. 2d) 



where 



2 p" 



-I 



8x - J x 2 - a 2 (N 2 - 1) 



6x + | x 2 - a 2 (N 2 - 1) 



> e dx , 



(4. 3a) 





2 (5)= [ 





2 2 2 


6x ^ 


x - or (N -1) 



► e dx , 



Sx + I x 2 - of 2 (N 2 - 1) 



(4, 3b) 



a = 2hJ3 = 2h(2ir/\) , N = € - is , s = 60 \ tr(mhos /m) , 



(4.4) 



and h denotes the antenna height in meters above the ground,, 

The variable x in (4. 3) assumes purely imaginary values from the lower limit to zero, and 
from there ranges over real values to the upper limit. The sign of the square root is chosen so that 
the real part is greater than zero. A derivation of the dipole impedances together with asymptotic 
expressions of (4. 2) is contained in the appendix. 

Ground pr oximity loss in terms of the input impedance as given by (4. 1) may be expressed as 



\ r = 10 log(r/r f ) = 10 log 



1 + Re(AZ/r f ) | , 



«4 e 5) 
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where Re (AZ/r^) denotes the real part of (AZ/r^), Using the relationships of (4. 2), ^ is 

plotted versus 

s = 60 X (r(mhos/m) 

in figures A-l to A-28 for various value s of and a . Notice that for a greater than about 

5, i. e„ , h/X ^ 0.4, the ground loss is always negligible; however, as the antenna height in wave- 
lengths approaches zero, L may become quite large for certain combinations of ground con- 

t, r 

ductivity and frequency. It is often possible to reduce the effect of the ground through the use of 
appropriate ground screens. Discussions of this subject are contained in papers by Wait and 
Surtees [ 1954] and by Wait [ 1956 ] . 

To show the variation in the normalized input impedance change for both the resistance 
(Re AZ./r^) and the reactance (ImAZ/r^), equations (4. 2) are evaluated and plotted in figures 
A-29 to A-64 as functions of | N | , cj>/2, and a where, 

/ 2 z\ ^ i 
|N| = ^ + s ) , cf>/2 = (%) tan" (s/^) , (4.6) 



with a and s as defined in (4. 4). It s hould be emphasized that figures A-29 to A-64 do not 
show the total input impedance, but only the variation due to the presence of the ground,, Impedance 
curves such as those in the figures are useful not only in estimating antenna requirements for lunar 
communication systems , but also may be used in determining the ele ctr omagnetic characteristics of 
material composing the moon. 

This section will conclude with a brief discussion of the Beverage wave antenna f Beverage, 
Rice, and Kellogg, 1923 ] . If considered for semi -permanent installations on the moon, it has the 
advantage of relatively s imply construction and operates most efficiently over poorly conducting soils. 
The wave antenna in its simplest form consists of a long horizontal wire situated a short height above 
the ground and terminated at one end through its characteristic impedance. It is a uni-directional 
antenna with the maximum gain in the direction of the antenna axis and toward the terminated end. 

Wait [ 1954 ] has shown that the vertical electric field component of a wave antenna is pro- 
portional to a complex factor T , termed the "wave tilt" and defined by 



T 

v 



(€ - 1) 

r 



/U-is) , (4.7) 



r 



and a function S* which depends on the wavelength X, antenna length i, an angle <j>' measuring 
the direction in which the antenna is pointing, and the propagation constant T defined as 

r ss a' + i|3m , (4. 8) 

where a % and m are real and (3 = 2tt/X « 
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For low antenna heights h an approximate expression for the values of a and m has 
been given by Carson [1926 ] : 



jr'/p c± m - 1 •. 



2 3/2 {(3h) N Jc r 2 +s 2 in(4h/d a ) 



(4.9) 



where d is the diameter of the antenna wire For an uninsulated infinite s imally thin antenna 

a 

actually lying on the ground, F may be expressed as [ Coleman, 1950 ] 

— — ^ 

(€ + 1) + S + ( € +1) 

\ r / r 



F s a 1 + ipm = (3 J m 2 - (^)(e + 1) + i({3/2) 



(4.10) 



The function S* , obtained by integrating the contributions of all the elements along the 
antenna, is given by 

q» - 1 - exp [ - (r - i(3cos , 

b _ — — (r _ ipcoSc p- }i - i^ii) 

It can be seen from (4. 10) that the propagation constant of a wave antenna lying on the ground 

approaches its free space value i (3 if the relative dielectric constant is very near unity and the 

conductivity is extremely small; in this case S 1 approaches unity a 

The wave antenna power gain [ Martin and Wickizer , 1949 ] referred to an isotropic antenna, 

(p /p. ), may now be expressed as 
wa iso 

(p /p. ) = ((3icos^) 2 | S' | 2 o | T | 2 (120 R /| Z | 2 ) , (4. 12) 

wa iso ' v 1 o 1 o 1 

where Z = R + i X is the characteristic impedance of the wave antenna. The factor I T I , 
o o o 1 v 1 

which was discussed in section 3 and defined by (3. 2), is shown in figure 2, and | S* | is graphed 
in figure 15 as a function of the two parameters (ai) and ((31) (m- cos 4> f )„ An example of the wave 
antenna pattern, 

|S| = (picos <\><) • |S' | , 

as a function of the direction angle 4> ? is shown in figure 16 for several values of (31 = 2ir(i/\) and 
for m = 1, a - 0. It is seen that the antenna is highly directional in the forward direction (cj)' = 0); 
other values of m and a modify the pattern somewhat, but the general shape is the same. 

A modification of the wave antenna which is also discussed in the paper by Wait [ 1954 ] is 
one in which the horizontal wire of length I = 2 b is center -fed and the ends of the antenna are not 
terminated with a load impedance. Assuming a current distribution along the wire of 

I(y) = I o [sinh(rb)]" 1 sinh jr(b- |y| )| , (4.13) 
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where I and y are respectively the current at and the distance from the feed points an integra- 
tion of the elements over the length of the antenna yields for the antenna pattern | S j : 



|S| = Ipbcos 4>' | • |S '| 



z\ coshx - cos (x_ -x_) T 2 I 

L 1 2 3 J L 



X l + V X 2 " X 3 



2i coshx ~ cos (x +x_) 
X L * " X 



2 / \ 

"i + ^2 + x 3 y 



2 2 
X l + X 2 



cosh cos x^ 



{4, 14) 



where x^ = a b, x^ = (3 b m , and x^ = (3b cos tj> „ 

Graphs of (4„ 14) are shown in figure 17 for the case of a = 0, which corresponds to very 
poorly conducting ground. The value s of | S | shown on the figure, however, should be used only 
with antenna lengths I such that the quantity (a I) is very near zero; if this is not the case, j S J 
should be calculated from (4. 14). The center -fed wave antenna power gain corresponding to (4„ 12) 
is given by 



p /p. 
wa iso 



= S 



120R / Z 
o 1 o 



i 2 ) 



(4.15) 



where T , R , and Z have been defined previously,, 
1 v 1 o o 

The theoretical de te rmination of L or L for a wave antenna in terms of its input 

t r r 

resistance is quite difficult. However, considering these quantities ess entially as power ratios 
measuring the effectiveness of the antenna, an estimate of their magnitude may be obtained from 
(4. 12) by assuming that 



L - G + 
t, r t, r 



10 log(p /p. ) . 
° \ wa iso / 



(4.16) 



5. NOISE EFFECTS 

The operating noise factor f [ Norton, 1962 ] appearing in (2, 4) provides a measure of 
the effects of noise on the receiving system (in previous papers [ Vogler, 1963 a, 1963b ] f is 
called the effective noise figure). It is defined by 

f=[P a /(k B t Q b)]/[r(s/n)] , (5.1) 

where p^ is the signal power available from an equivalent loss -free receiving antenna, (k^ t Q b) 
is a reference Johns on -noise power [ see (2. 3) ] , and r(s /n) is the receiver predetection output 
s ignal -to -nois e ratio. 

Since p is the power available from a loss -free antenna, its relationship to the power 
available from the actual antenna, p^ , may be expressed as 

P : , =i P„ . (5.2) 

a c r 
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where 4^ is the antenna circuit loss factor. Because the present study is primarily concerned with 
low antennas over poorly conducting ground, the assumption is made throughout that 4 £h 4 , 
where is the ground proximity loss discussed in section 4 

In terms of the noise factors of the receiving system's component parts, f can be expressed 
in the following form: 

f = f e- 1+f c + i c<V )) + 'c V^- 1 ' • < 5 - 3 > 

where f^, f^, f , and denote respectively the noise factors of the external noise, antenna 

circuit, transmission line, and receiver,, In terms of the loss factors 4 and 4 „ of the antenna 

c ti 

circuit and transmission line (i. e , the ratio of the available input to output powers of the component), 
and f are given by 

f c = 1+ ( V l)(t c /t o ), £ u =H(i a -l)Vy . (5.4) 

where t^ and t^ designate the absolute temperatures of the corresponding components and t is 
a reference temperature. 

Generally speaking, the factors f and f will predominate in calculating values of f , 
especially at the lower frequencies [ Crichlow, Smith, Morton, and Corliss , 1955]; at high fre- 
quencies f depends more on f and f , which are best obtained by direct measurement. By 
arbitrarily assuming the transmission line loss to be negligible (i = 1 ) and the antenna circuit 
temperature very near the reference temperature, ( 5. 1 ) may be written as 

f = f + 4 f - 1 , 4=1, t = t . (5 5) 

e c r t4 co 

With the appr oximation 4 £^ 4 , this then becomes 

c r 

£ ca f + 4 f - 1 , (5. 6) 

err 

where 4 has been defined in section 4. 
r 

The principal sources of external noise received by an antenna located on the moon are 
( 1 ) galactic or cosmic, (2) earth-based, and (3) solar „ The amount of noise af f e cting the signal 
depends on the location, orientation, and directional gain of the particular antenna under considera- 
tion,, Some preliminary studies [ Page, 1962 ] of earth -based and solar noise components indicate 
that, in the case of a quiescent sun, solar noise is negligible compared with galactic noise unless 
the antenna beam is oriented directly toward the sun, On the other hand terrestrial atmospheric 
noise might become comparable to the galactic noise in the lower frequency range; however, this 
assumes a negligible ionospheric shielding effect 

Galactic noise levels have been fairly well established at the higher frequencies through 

earth-based measurements . Figure 1 8 shows an estimate of the external noise factor F = 10 log f 

B e s e 

obtained from galactic noise data compiled by Menzel [ 1961 ] and by Hartz [ 1963 ]; the reference 
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temperature assumed here is t Q = 288. 37 K. Empirical expressions used to plot the galactic noise 
curve are 



f = 6o467 X 10 6 f " 3 , f > 200 , 
e mc mc 



f = 1 585 X 10 5 f " 2,3 , 10 < f < 200 , 

e mc mc 



(5. 7a) 



from the Menzel data, and 



4 - l.£ 
f = 5.012 X 10 f 

e mc 



0.5 < f < 10 

mc 



(5.7b) 



from the Hartz data obtained from Topside Sounder measurements, The upper and lower dashed lines 

2 3 18 

on the figure are extrapolations of the f " and f ° frequency variations. 
& r mc mc 

It is physically apparent that at some lower frequency the galactic noise curve will at least 
level off. Until information from low frequency measurements are available, the extrapolations of 
figure 1 8 may be used to provide (hopefully pessimistic) estimates of galactic noise. 

6. POSSIBLE IONOSPHERIC EFFECTS 

At the present time there is no indication that the moon possesses a neutral atmosphere of 
any significant magnitude Because of this , it is unlikely that lunar surface propagation will be 
affected by ionospheric considerations to the same extent as terrestrial propagation. However, 
studies have been made showing the possible existence of a weak lunar ionosphere caused by solar 
wind. This section will discuss briefly the effect on propagation of two forms of possible ionospheres. 

If a negligible magnetic field is assumed and the motion of heavy ions is neglected, the 
ionosphere may be characterized by its index of refraction KL as follows : 



N 2 = 

1 



1 « 



r 2 2 
f + f 

mc v 



f f 
v cr 



( f 2 + f 2 



mc \ mc 



(6.1) 



where f denotes the signal frequency, f the collision frequency, and f the critical or 
mc & v cr 

plasma frequency, all measured in Mc/s. In terms of the electron density/ cm , n, the critical 
frequency is given by 

f Z m 8JX 10" 5 n, (6.2) 

cr 

where, in general, n is a function of the height h above the ground: n = n (h) . 

One form of ionosphere that may possibly occur extends from the moon's surface to some 

distance above, the electron density decreasing with height. Elsmore [ 1957 ] , using studies of 

lunar occultation of the Crab Nebula, has hypothesized surface electron densities on the order of 
3 4 3 

10 to 1 electrons /cm , For signal frequencies near or less than the collision frequency, this 
ionospheric model presents a propagation problem which, as far as the author is aware, has never 
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been investigated and is beyond the scope of the present paper. Considerable theoretical study- 
would be necessary to determine propagation characteristics at these low frequencies. 

For radio frequencies such that f /f < < 1, an estimate of the effect of a surface ion- 
^ v mc 

osphere on propagation at frequencies greater than the plasma frequency may be obtained through 
the concept of effective radius [Burrows and Atwood, 1949 ] . The effective radius, r^ , intro- 
duces a correction to the actual radius, r , that allows for the refraction of the radio wave as it 

o 

travels through the medium (in this case the model ionosphere) above the ground. Expressed in 
terms of an effective radius factor k [ Wait, 1961 ] : 



1 + 



r N. (h=Q) 
o i 

N (h=0) 



(6.3) 



where N. (h = 0) denotes the gradient of the refractive index with respect to height evaluated at the 

surface o 

With the restriction that f /f < < 1, (6. 1) may be rewritten as 

v mc 



N 2 « 1 - f 2 /f 2 = 1 - 8,1 X 10" 5 n(h)/f 2 
l cr mc mc 



(6.4) 



and k then becomes 



1 + 



4.05X 10 r {-n'(0) } 

o 

N 2 (0) - f 2 
l mc 



1 + 



0.07 { -n (0) } 



f 2 - 8.1X 10" 5 n(0) 
mc 



.1 



(6.5) 



where n(0) is the electron density /cm at the surface and n (0) is the electron density gradient 
with respect to height in km evaluated at the surface. It can be seen that (6. 5) may be used only 
for radio frequencies such that 



-3 



f > f = 9 X 10* N/n(0) ; 
mc cr 



thus, for an electron density at the surface of, say, n(0) =10 , the effective radius concept should 

be restricted to frequencies of 1 Mc/s or greater. 

The curves of figures 9-14 will now provide estimates of the attenuation A if they are 

read as functions of a modified x ' , x '(k) , and a modified K. , K (k) , where 

o o h, v h, v 



x (k) = k 3 f 3 d (km) , 
o mc o 



K(k) = 3.02 X 10" 2 k"y | T 



Q 



(6.6) 



and k is given by (6. 5). Notice that x (1) and K(l) are simply those parameters discussed in 
section 3 and defined by (3. 5) and (3. 1); the effective radius factor k does not modify b^ ^ 
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To show the variation of k with frequency, (6, 5) is plotted in figure 19 for several values 
of the surface electron density gradient n'(0). The values chosen do not represent any experimental 
evidence but are used merely to show the variation of k with the gradient. Figure 19 shows that, 
as the gradient becomes more negative, k decreases and, in general, a higher attenuation of the 
signal is to be expected,, Should a positive gradient exist, less attenuation of the signal would occur. 

Another possible lunar ionosphere that has recently been investigated [Weil and Barasch, 

1963 ] is similar in form to the terrestrial ionosphere Weil and Barasch, through theoretical 

considerations of solar wind influence near the moon together with solar stream data, postulate an 

3 

ionized elevated layer with an electron density of about 350/cm at o 6 lunar radii above the 
surface . The layer occurs only in the general direction of the sun, and their numerical values are 
based on quiet sun data c 

The theory of this form of ionospheric propagation is well established and, given the neces - 
sary data, will provide good estimates of propagation capabilities. The remainder of the section 
will discuss briefly the simple geometrical -optics approximation to ionospheric propagation,, A 
detailed account of the subject is contained in a paper by Wait [ 1962 c ] . 

In the geometrical-optics approximation, the field at a point is considered to be the sum of 
multiply -reflected radio rays arriving from the transmitting source, the reflections occurring 
between the ground surface and the ionosphere. At distances less than its caustic, the single hop 
component of the field, i. e. , the ray undergoing only one reflection from the ionosphere, will be 
the most significant. However, the total field is given by the sum of all the rays added together in 
proper phase; at short distances from the transmitting antenna, the ground wave component must 
also be included. 

In order to gain some idea of the effect on lunar propagation of an elevated ionosphere, the 
attenuation of the single -hop sky wave is evaluated, assuming the Weil and Barasch model iono- 
sphere with a sharp lower boundary. The assumption of a sharp boundary rather than a more 
realistic one will give a lower estimate of the attenuation than is actually obtainable , but will still 
provide a reasonable estimate for the purposes of the present study. Letting 0{ = d fx ) , 
p( = h/ r^), cj>^ f and (j>, denote the geometrical quant it ie s indicated in figure 20, the sky wave 
attenuation A g may be expressed as 



A ~ -20 log 
s to 



C sin 2 § \ z 

2 0,/9 g /' l 1+R gl ° l R il 



1' 



(6.7) 



where R^ and R„ are Fresnel reflection coefficients of the ground and ionosphere respectively : 



N 2 cos 6 - J N 2 - sin 2 cb 
g g g r g z 

R = _ — — n = e - i s , (6.8) 

g z \ 1 ^ g g g 

N cos cb + J N - sin cb 
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Figure 20, Geometry for single -hop sky wave propagation. 
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2 2 2 

N. cos 6. - \ N. = sin 6. 

1 T l N 1 '1 



R. = 



2 2 2 

N. cos cb. + IN. - sin <J>. 

i i \j i T i 



N. =€. - is. , 

ii i 



(6.9) 



and the angles cf> 



and c|). may be calculated from 



c|> = tan 



(1+p) sin(9/2) 
(1+p) cos (9/2) - 1 



j>. = sin 
i 



i+p 



(6. 10) 



The quantity C in (6, 7) is a convergence coefficient defined by 



C = 



1+p 
" cos (9/2) 



1 + 



sin (9/2) 
(1+p) cos (9/2) - 1 



(6. 11) 



and 9 ^ is related to the total distance traversed by the ray: 



9 = Jp 2 + 4(l+p) sin 2 (9/4) 



(6,12) 



Notice that (6. 1 1) [ and therefore (6. 7) ] is valid only as long as 



< 9 c = 2 cos" [ l/(l+p) ] 



9 corresponds to the distance at which a caustic is encountered, and other formulas must be used 
c 

to calculate the attenuation in this region,, The distances considered here will be only those such 

that 9 < < 9 , 
c 

At VLF and for the distances of interest, the reflection coefficient of the ground R^ is near 
unity even for the very poorly conducting ground expected on the moon, On the other hand, the ion- 
ospheric reflection coefficient R can have a pronounced effect on propagation estimates . Figure 
2 1 shows a plot of | R^ | as a function of the ionosphere height p = h/r^ and 



s. = f f "/f ( f 2 + f 2 
i v cr mc \ mc v 



obtained from (6. 1)„ Notice that | R^ | depends quite critically on the parameter s^ and, hence, 
on estimates of n and f^„ If an elevated lunar ionosphere does exist, data concerning these quan- 
tities are necessary for VLF propagation predictions. The purely geometric factor 
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C sin 2 cj>_ )/{ 2 8 1 /9 

appearing in (6. 7) depends only on the distance d^ between the transmitting and receiving points 
and on the height of the ionosphere; curves of this factor are presented in figure 22 for various 
values of p • 

The attenuation A g of a single -hop sky wave relative to a free space inverse distance field 

is shown in figure 23 as calculated from (6. 7) a For these curves the ground is described by 

2 2 
N = 2 - i 180, and the ionosphere at a height p = 0.6 is described by N. = 1 - i s. , where s. 

g r ill 

assumes the value s 1 and 10„ At a radio frequency of f = 0.01, for example, the curves 

. mc 

—4 

would be valid for a ground conductivity = 10 mhos /m and an ionosphere described by two 

values of the ratio (n/f ) cs* 125 and (n/f ) as 1250 ; the latter values assume f >>f , The 

x v v v mc 

dashed curves in the figure show the vertically polarized ground wave attenuation for 

f - 0.01 and 0. 1: this component of the field is dominant at the shorter distances. It can be 
mc 

seen that over the range of distance shown, the difference in attenuation between the two sky wave 
curves is about 10 db 9 this difference arising from the ionospheric reflection coefficient R . 
A change in ground constants at VLF produces very little difference in attenuation for a given s^„ 
According to Weil and Barasch [ 1963 ] the ionosphere would exist (with any appreciable 

o 

magnitude ) only within about 60 from the sun-moon line; thus, it would not be a dependable 
mechanism for long-term propagation. It should be remembered, however, that it could occasion 
interference effects in VLF ground wave c ommuni c ati on circuits at those distances where the ground 
wave and sky wave components are comparable in magnitude. 

7. CALCULATION OF REQUIRED POWER 
Attempts have been made recently to deduce the electromagnetic properties of the moon's 
surface material through the use of radar data [Senior and Siegel, I960 ] . Although differences 
exist concerning the exact interpretation of the data [ Brown, I960; Daniels , 1961 ] , there is 
general agreement that the relative permittivity is not far above unity and that the conductivity is 
quite low. In any cat'e, most of the graphs discussed in the preceding sections are applicable to a 

wide range of values of e and <r « 

r 

Senior and Siegel estimate the relative dielectric constant and conductivity of lunar surface 
material to be = 1. 1 and cr = 3.4 X 10 mhos /m. Materials constituting the earth's crust 
have been found with conductivities of this order, but not with such a low permittivity. For the 
purpose of arriving at some estimate of lunar propagation conditions , it will be assumed that 

ranges from 1.1 to 2.0 and g- lies between 10 ^ and 10"^ mhos /m . 

o 

By choosing a reference temperature t = 288. 37 K (as in section 5), the required power 
given by (2. 4) may be rewritten as 

P = L - (G + G ) + L + R + F + B - 204 (7. 1) 

t p t r* t 

where the propagation loss L^ for ground wave attenuation between isotropic antennas is given by 
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L = 32,45 + 20 log d (km) + 20 log f 



mc 




(7.2) 



For sky wave propagation, JL would be calculated by (7. 2) after replacing A by the A dis- 

p t s 

cussed in section 6. 

Figure 24 shows curves of for a vertically polarized ground wave over homogeneous 

ground as a function of frequency and distance. Values of- A were obtained from figures 11 
through 14 by linear interpolation in- both the and b^ directions, being given by figure 2 

and (3, 1) with | Q | = 1, and b being read from figure 4. Notice that varies inversely 

with conductivity at the lower frequencies, while at high frequencies the effect of cr variation 
becomes negligible. 

To show the use of (7. 1) in estimating required power for a particular communication 

system, the following example will be considered. Assume a communication system consisting of 

a wave antenna lying on the moon's surface and transmitting towards a short vertical electric dipole 

t ° 

placed some distance away and in an optimum direction from the transmitter (§ ^ ). An estimate 
of L - G may be obtained from (4. 16) and (4. 12) by assuming the free space value for the 
wave antenna characteristic impedance: 



For the case of an antenna of length I = \/4, L - is plotted versus frequency in figure 25, 
Again, the antenna loss varies noticeably with conductivity at the lower frequencies, but not at the 
higher ones. 

A comparison of figure 18 and figures A- 1 and A-2 shows that for a dipole at a height of, 
say, h = X/ 16 the external noise tends to "blanket out" the effect of the antenna circuit loss,, This 
is shown in figure 26, which is a plot of the operating noise factor as given by (5. 6) with an assumed 
receiver noise factor of = 4. Now by using figures 24, 25, and 26, and the dipole gain of 

= 1.76, the required power delivered to the transmitting antenna terminals for a given type of 
service may be estimated from (7. 1). Figure 27 shows values of P - (R+B) as a function of 
frequency for antenna separation distances of 10, 100, and 500 km . 

To estimate the required power delivered to the terminals of the transmitting antenna, the 
type of communication service desired must be designated, thus specifying R and B, For 
example, with standard broadcast service and a bandwidth of 10 kc/s, R may have the value 
39 db [CCIR, 1951 ] and B = 40 db so that 79 db should be added to the curves of figure 27 to 
obtain the required power in decibels above 1 watt. Thus the power required for this type of service 
in the example system at a range of 10 km and for a frequency of 300 kc/s would be about 10 db 
or 10 w. For a low-grade voice communication service and 6 kc/s bandwidth, R = 9 db [ CCIR, 
1951 ] , B = 38 db, and 47 db would be added to the curves. The required power in this case at a 
distance of 100 km and for a frequency of 100 kc/s would be about 16 w- a 

If fixed wave antenna lengths and dipole heights are considered, rather than the variable 
lengths and heights assumed for figure 27, ground proximity loss curves may be recalculated from 
the equations and graphs of section 4 and required power estimates obtained from (7. 1). Figure 28 
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Figure 27. P -(R+-B) for wave antenna and vertical dipole over 
homogeneous ground. 
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shows a plot of the quantity L - G for fixed wave antenna lengths of 1 = 100 m and 1000 m 

t t -3 -4 

and ground conductivities of or = 10 and 10 mhos/m with = 2. It can be seen that, at 
the lower frequencies, the loss associated with the 1000 m antenna is about the same as that shown 
in figure 25; however, a higher loss occurs at the lower frequencies if the 100 m length wave 
antenna is used. 

The ground proximity loss for a short vertical electric dipole at fixed heights h = 10 m 
and 100 m is shown in figure 29 for the same ground constants assumed in figure 28. As before, 
the extrapolated value s of the external noise factor of figure 18 will tend to override the dipole loss 
at the lower frequencies, but it should be emphasized that this loss will contribute significantly to 
the calculation of f at low frequencies if future measurements show f to be considerably lower 
than the extrapolated values assumed here 

Using figures 24, 26, and 28, the quantity P - (R + B) is shown in figure 30 for the case of 
a wave antenna of length I = 100 m transmitting in the optimum direction towards a short vertical 
electric dipole at a height h = 10 m e For the example of standard broadcast service used pre- 
viously, the required power P^ at a range of 10 km and for a frequency of 300 kc/s would be 
about 12 db or 16 w„ For the low-grade voice c ommunic at ion service the required power at a 
distance of 100 km and for a frequency of 100 kc/s would be about 30 w. At high frequencies 
the estimates of P shown in figure 30 will be somewhat pessimistic, since a small height gain 
effect will occur due to the height of the dipole. Graphs for calculating this effect may be obtained 
from another paper [ Vogler, 1964]. 

The combination of a wave antenna and dipole and the particular lengths and heights used in 
the preceding example were arbitrarily chosen. As pointed out in section 1, the purpose of the 
c ommunic ation system together with the necessary restrictions placed on its physical components 
must first be established before specific power requirements may be estimated. Combinations 
of antennas other than the one assumed in the present example should be investigated before a 
decision is made as to what will constitute an efficient lunar communication system. 

In the particular idealized system described by figures 27 and 30, a number of points should 
be noted: ( 1) the required power does not vary appreciably over the range of assumed, but 

does depend significantly on the conductivity at low frequencies; (2) propagation out to s ome what 
beyond 100 km is practical for most types of systems and service, at least at MF or below; 
(3) the curve s indicate an optimum frequency exists, depending on the conductivity of the lunar 
surface and the range of propagation; for the model assumed and for the c onductivitie s and dis - 
tances shown, the optimum frequency lies in the LF band. It should be kept in mind that the oper- 
ating noise factor of the receiving system was considered to be a function only of the galactic noise 
(which, of course, is extrapolated at the lower frequencies) and a rather low receiver noise factor. 
If the receiving antenna loss i were large enough, it is apparent from (5. 3) that F would have 
higher value s than those assumed. Also, during the lunar day, the antenna temperature t^ 
would increase, thus making the operating noise factor even higher. 
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Figure 30. P t «fR+Bj for fixed length wave antenna and fixed 
height VED over homogeneous ground. 
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8, CONCLUSIONS 

A study has been presented concerning various aspects of lunar surface communication; it 
is hoped that these results may serve as a basis for future investigations of the subject, An attempt 
has been made to present the main curves in a generalized form such that, as changes and refine- 
ments in our knowledge of the moon arise, the curves may still be useful in predicting lunar 
communication requirements. 

It is apparent that further studies, both theoretical and experimental, would be useful in the 
following areas: 

1, A more precise evaluation of lunar ground conductivity, together with more information 
concerning layering (if any) of the materials composing the moon, 

2, Further investigation into the structure and magnitude of possible lunar ionospheres , 
especially the dete r mination as to whether an ionized layer exists at the surface, and if 
so, its effects on propagation at the lower frequencies, 

3, Further studies of noise effects including terrestrial noise sources, noise during high 
solar activity, and low frequency galactic noise, 

4, Investigations of various antenna combinations , taking into account antenna ground losses 
for those situations in which a ground screen will be impractical, 

5 e A study of knife-edge diffraction and obstacle gain effects caused by irregularities in the 
moon's surface structure, 
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10. APPENDIX: INPUT IMPEDANCE OF DIP OLE ANTENNAS 



It is assumed that an elementary dipole of length dl is located at the origin of a cylindrical 
coordinate system (p , 9, z) at a height h above an isotropic and homogeneous conducting half- 
space characterized-by a dielectric constant and conductivity g-y The interface between 
the two media is the plane z =. - h 9 and for z > ~h the dielectric constant is assumed to be equal 
to the free space value The magnetic permeability of the whole space is taken as the free space 
value jj, (MKS units are used here),, In the case of a vertical dipole, the antenna is oriented in 
the z direction, whereas for a horizontal dipole the orientation is in the 8=0 directions The 
elementary magnetic dipole may be considered equivalent to a small circular loop with its axis in the 
direction of the dipole and with an area dA = &l/p^» where (3^ is the free space wave number 
( = 2tt/X.) For practical short dipoles having a linear distribution of current, 61 should be replaced 
by an effective length = di/2 

The ele ctr omagnetic field of an elementary dipole above a conducting plane may be obtained 

from the electric and magnetic Hertzian potentials, n and the components of which, if the 

vectors are referred to a cartesian coordinate system (x, y, z), satisfy the scalar wave equation 
2 2 

( V + (3 ) 3TI = . Using the cylindrical coordinate form of the wave equation, particular s olutions 
are found by the method of separation of variables resulting in expressions for the IT components 
involving combinations of the eigeiif unction 



where X denotes the eigenvalue (following Sommerfeld's notation), J is a Bessel function of the 
first kind and order n, and the C ' s are functions determined from the boundary conditions . These 
conditions in the present instance require that the tangential components of the electric vector E 
and the magnetic vector B be continuous at the interface between the two media (air and ground), 

In the case of vertical dipoles, E and B may be expressed in terms of Hertzian potentials 
having a z -component only Thus in the uppe r medium 



where the superscripts d and r refer to the direct or primary stimulation and the reflected or 
secondary stimulation, respectively. After solving for the coefficients C through the boundary 
conditions, one finds that for 




(A-l) 




(A-2) 




o 



o 



(A-3) 



n r = 

ez 




e 



-(z + 2h)u 

o 



(\/u )dX, 



o 



49 



VMD: TI 



iK 

, d m 



mz p 



J o (X,p) e (\/u o )d\, 



(A-4) 



iK 



n 



mz (3 



n r u~u.. -(z + 2h) u 



u + u 1 

o 1 



where 



(A-5) 



P q and (3^ are propagation constants of the upper and lower media, respectively, and the K's are 
amplitude factors given by 



P I (di) 

K = -2-2 , K 

e 4ttco€ m 



3 jx I (dA) 
o o o 

______ 



(A-6) 



For horizontal dipoles, the Hertzian potentials require both an x- and a z- component. Thus, 
corresponding to (A-2) , 



n, = n x + n z = m +n ) x + n z 

hor x z xx z 



(A-7) 



After solving for the C's through the boundary conditions, one obtains : 

iK 



HED: n 



r o o 



• T jr OO 

iK n - u - u 

n r = f J (Xp) L I e 

ex (3 J o x Vl 1 u + u J 
o o v. o 1 J 



(z + 2h)u 



(X/u )dX, 



iK cosG r 2(u - uj ^ -(z + 2h)u _ 

_|_ C J( x p > J_^L_L L e V/^dX, 

P o J o 1 Ln u +u, J ° 

o 1 



(A-8) 



HMD: n 



i K «°° - z u 

= — — \ J (Xp) e °(\/u )dX, 
mx 6 J o o 
r o o 



K r N u - u. , -(z + 2h)u 

m f J o (Xp) j , ° e °(X/uJdX, 



mx (3 



o o 



iK cos ,-,<» 
m 



N u + u. 

o 1 

2{u 



J J ( xp) {-2-2-1} 

^O U + U„ J 



(u ~u, k ~(z + 2h)u 7 

° 11 °{x 2 /p o 2 )dx, 



o 1 



(A^9) 



where u , u,, N , K , and K are defined in (A-5) and (A-6), 
o 1 e m 
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With the Hertz vector II determined, the electromagnetic field in the upper medium 
derived through the following relationships : 



E = grad div II + p n ^ , (electric dipoles) 



B = grad div n + p II . (magnetic dipoles) 



The components (in cylindrical coordinates) of interest in this paper are found to be 



,2,2 2 \ f Z , Z Z V d r 

VED: E (p, z) = [d dz + p II = [d /dz + (3 )(n + n 

~ r v r o / ez \ r o / \ ez ez 



iK 



(3 J o 
r o o 



J (X.p) 



z u r Nu - u ^ -(z + 2h) u 



e ° + 



N 



— } 

u +u, J 



(X /u Q )d\, 



o 1 



HED: E (p,z) = (9/9p) div n + (3 cos 9 n 

p e o ex 



iK cos 

e 



p r / 2 2 f J l l ^ } 1 ^ 

J o [(e j o< x '>- x i j o< x p)-— }; 



r J (\pk n f u -u -(z + 2h)u 

+ 1 P^J (Xp)-xNj (\p). 1 



\p f / I u +u 
r J o 1 



2 2. 



+ U Up) r—Ji —r- 2 f e 

V ° Xp y l p 2 (N 2 u +VL) J 

o o 1 



(\/u )d\, 



/ 2 2 2\ f Z , 2 2 N \/"d 

VMD: B (p,z)= a /3z +8 ) n =(d dz + p ) i n + n 
z r V r o / mz V r o / V mz 



r 

mz 



iK p oo 



. -zu ,- u - u - -(z + 2h) u 
o f o 1 J 



° 1 (\ 3 /u Q )d\ $ 



* 2 

HMD: B (p, z) = (d/dp) div n + S cos n 
p r r m o i 



iK cos n 

_5L c 



B Z J (\p)-\ 2 ^| J (Xp) - 



J-Up)^ »zu o 



+ (p o j o (x P ,-x |j o 



Up) 



J^P),\jN u ^ u n -(z + 2h)u Q 

_____ | _ 
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+ U <M- Xp 



JJXpK 2\ u (u -uk -(z + 2h)u 
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1 (\/u )d\« 

! o 



To calculate the dipole input impedance using the emf method [ Wait, 1953 s 1962 d; King, 
1956 (p. 258) ] requires a knowledge of the axial field at the antenna. This may be obtained by taking 
the limits of (A- 11) as p and z approach zero (and with 9 = in the case of horizontal 
dipoles). The impedance Z is then given by 



VED: Z = Lim4E (p, z) di/l 

p-^o 1 z 
z->o 



iK di 



00 ^o l. ^ N u + u, 



N u - u, -2h u 



o 1 



(X /u )dX, 



(A-12a) 



HED: Z = Lim Je (p, z) dl/l L 
0—0 1 P 



p— o L p 

z— *-o 



iK di 
e 



(3 I 
o o o 



I 



^ o 1 



2 2 

r X u (u - u, ) s -2hu . 

I o o 1 I o 

+ j — — - — — - y e 

L 6 (N u + uj J 

r o o 1 



(X/u )dX, 

o 



(A-12b) 



VMD: Z = Lim \ coB (p, z)dA/l 

p— o 
z-*o 



icoK dA ° 
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(3 I J 

r o o ^ 



L l U o +u l 



u - u . -2hu 



(X /u )dX, 
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(A-12c) 



HMD: Z = Lim J co B (p, z)dA/l 

p— o LP 



i co K dA rt oo 

m r 



^ N u + u 

o 1 



2hi 



2 2 
X u (u - u ) 



■ (3JN u +1^) J 



-2h\ 



(X/u ) dX« 



(A-12d) 
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Notice that the effect of the ground is contained in terms involving the factor exp (- 2hu ); thus, the 

o 

change in impedance AZ due to the presence of the ground may be defined as 



AZ = Z - Lim Z = Z - Z f 

h-*oo 



(A-13) 



and the free space antenna resistance r is obtained from the relationship: 



r = Re J Lim Z 



(A-14) 



With the use of (A-13) and (A-14), one can now derive the expressions for the normalized 
input impedance change AZ / r f given by equations (4c 2) of section 4; here, the change of variable 



f i ? 
x = 2hu = 2h J \ - 6 

O N r o 



has been made and the following relationships used: 



(\d\/u Q ) = (dx/2h), a = 2hP Q , u x = ^lx 2 -« 2 (N 2 -l)/2h, \ 2 = ( a 2 +x 2 )/(2h) 2 . (A-15) 

Limiting expressions of (4« 2), derived from integrations by parts of the integrals involved, 
are given by the following e quations For orN < < 1 : 



.f Fi(0)+ f^lY F (D) ] +i (iil±|^) e- A , 

L J ^ 2n> J 



HED: AZ/R r c~ 



3(N + 1) 



f "~ 8 



(A- 16 a) 



(A- 16 b) 



VMD: AZ/R f c 3(N / ^ 



Fi (0, + (N+1; F 2 (0) 



(A-I6c) 



HMD: AZ/R,^ 



3(N + 1) 



"f " 8 



Fi (D)+ ( iij+l) F 2 (0) 



, . f 3(N + 5)D ] -A , A _ , 

+ 1 { 16 a j e ' (A - l6d > 



where 



F^X) = X + {1 - X(l + X) } d + |(1 + X)(l - X 2 ) d 2 , 



X 



F 2 (X) = - ~ + {1 + X(l-X)}d +{1-X(l + X)(2-X)}d 2 + i{l + X(l-X)(2-X 2 )}d 3 - I (1 - X) (1 - X^ d 4 , 

(A-18) 

iA-19) 
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N + 1 
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For a N > > 1 : 



(A-20a) 



HED: AZ/R f ~ -%{|^} 



| f (l - of 2 + 2/N) + ilo? + -4- U 2N " 3 



N 



(A-20b) 



VMD: AZ/R. 



3 f N-l 
3 1 N+l 



(A-ZOc) 



HMD: AZ/R 



3 r N-i 

■f ~ - 3 N+l 



ie" lQf (A-20d) 



For the case of dipoles above a perfectly conducting plane (N-*°°), equations (A-20) reduce to the 
well-known expressions : 



VED: AZ/R f — , 3 



{sin a - a cos a} + i {cos a + a sin a} , 



(A-21 a) 



HED: AZ/R f J [~ {(1 - a 2 ) sine* - a? cos a} + i {(1 - a? 2 ) cos a + or sina? } ] , (A-21b) 

2q? i- 



VMD: AZ/R f — - ( -~ y [~ {sin a? - a cos a?} + i {cos a + a sin a} "j , 



(A-21 c) 



HMD: AZ/R 



f V ? 3 
2 a? 



2 2 
{(1 - a? ) sin a? - a cos a?} + i {(1 - a ) cos a? + a? sin a?} I • (A-21 d) 
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12. GRAPHS OF GROUND PROXIMITY LOSS, 1^ ^ = 10 log(r/r f ) 

for 

Vertical Electric Dipoles (VED) 
Horizontal Electric Dipoles (HED) 
Vertical Magnetic Dipoles (VMD) 
Horizontal Magnetic Dipoles (HMD) 



s = 60 \ cr(mhos/m), a = (2h) (2tt/M , 

€ ; relative dielectric constant of ground 

r 

cr{mhos/m): conductivity of ground 

h: height in meters of antenna above ground 

X,: wavelength in meters 
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S = SO Act ( mhos/m) = 1.8 x I0 4 cr (mhos/m)/f 



Mc 



Figure A- 17<* Ground proximity loss L for VMD« € = 5 • 
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13, GRAPHS OF NORMALIZED INPUT IMPEDANCE CHANGE, AZ/r f 



for 



■Vertical Electric Dipoles (VED) 



Horizontal Electric Dipoles (HED) 



Vertical Magnetic Dipoles (VMD) 



Horizontal Magnetic Dipoles (HMD) 




4»/2 = O^tan^s/^) , 



+ s 



s = 60 \ cr( mhos /m), a - (2h) (2ttA) , 

c : relative dielectric constant of ground 
r 

o-jmhos/m): conductivity of ground 

h: height in meters of antenna above ground 

X: wavelength in meters 
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Figure A-47. Normalized input impedance change for VMD, 
| Nj = 1.05 . 
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Figure A-49. Normalized input impedance change for VMD, 
|N| = 2 . 
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Figure A-52 c Normalized input impedance change for VMD, j N j = 5 
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Figure A«56. Normalized input impedance change for HMD, 
|N| = 1.05 . 
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Figure A°»58„ hozn-siiz^d input impedance change for HMD, 
■ N " - I 
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